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We consider the discrete Laplace operator on Erdos-Renyi 

random graphs with iV vertices and edge probability p/N. We are 
interested in the limiting spectral properties of A^ N ' as N — > oo in 
the subcritical regime < p < 1 where no giant cluster emerges. We 
prove that in this limit the expectation value of the integrated density 
of states of A^' exhibits a Lifshitz-tail behavior at the lower spectral 
edge E = 0. 

1. Introduction. The last decades have seen a growing interest in spec- 
tral properties of linear operators defined on graphs, mostly of the adjacency 
matrix or the graph Laplacian [10, 12, 13, 29]. The aim is to see how prop- 
erties of the graphs are reflected in properties of the operators and vice 
versa. 

Spectral properties of random graphs, however, still remain to be uncov- 
ered to a large extent. The mostly recent works [1, 3, 9, 23, 34, 35] deal with 
random subgraphs of an infinite graph, such as obtained by a percolation 
model. Their results range from ergodic properties of the spectrum to the 
existence and regularity properties of the integrated density of states, as well 
as its asymptotic behavior near spectral edges. 

A different prototype of random graphs was introduced by Erdos and 
Renyi [15], see also [4, 16] for more recent accounts. Here, one is interested 
in a scaling limit N — > oo of an ensemble of graphs with TV labeled vertices 
and an TV-dependent probability measure. The problem is to get spectral 
information on the Laplacian or other matrices associated with the graph 



Received April 2005. 

1 Supported in part by the Deutsche Forschungsgemeinschaft (DFG) Grant Mu 1056/2- 

1. 

AMS 2000 subject classifications. Primary 15A52; secondary 05C50, 05C80. 
Key words and phrases. Random graphs, spectra of graphs, Laplace operator, eigen- 
value distribution. 

This is an electronic reprint of the original article published by the 
Institute of Mathematical Statistics in The Annals of Applied Probability, 
2006, Vol. 16, No. 1, 295-309. This reprint differs from the original in pagination 
and typographic detail. 



1 



2 



O. KHORUNZHIY. W. KIRSCH AND P. MULLER 



in this limit; see, for example, [2, 17, 20, 24]. As compared to the situation 
described in the previous paragraph, this one here shares more similarities 
to the spectral theory of large random matrices, which was originated by 
Wigner [36, 37]. In these studies the primary questions are related to the 
existence and explicit form of the mean eigenvalue distribution function of 
N x iV-random matrices in the limit N — > oo or, in other terms, of the 
limiting integrated density of states. 

In the present paper we study a problem that joins the two branches de- 
scribed. We consider the discrete Laplace operator (the graph Laplacian) on 
Erdos-Renyi random graphs and show in Theorem 2.5 that the asymptotic 
behavior of its limiting integrated density of states at the lower spectral 
edge is given by a Lifshitz tail with Lifshitz exponent 1/2. This means that 
the occurrence of eigenvalues right above the lower spectral edge is a large- 
deviation event. The dominant contribution to the Lifshitz tail is provided by 
the linear clusters of Erdos-Renyi random graphs. We refer to Remark 2.6(i) 
and (ii) for a brief outline of our proof. Such a strong probabilistic suppres- 
sion of eigenvalues near a spectral edge was first quantitatively described 
by I. M. Lifshitz in the physics literature to account for certain electronic 
properties in disordered materials [26, 27]. It is nowadays a well-understood 
phenomenon in the mathematical theory of random Schrodinger operators 
[8, 22, 25, 30, 33]. 

2. Model and result. Given a natural number N >2 and a positive real 
p £ ]0, N[, we consider Erdos-Renyi random graphs with N vertices and 
edge probability p/N. There are many interesting phenomena when allowing 
p to grow with N, see, for example, [4, 15], but in this paper we consider 
the sparse case where p is fixed and does not depend on N. The graph Q( N ) 
is a random subgraph of the complete graph KS N ^ with N labeled vertices. 
Edges are distributed independently in with the same probability p/N. 
In other words, if i s anv given subgraph of K.( N ^ with M edges, then 
it is realized by Q( N > with probability 

(2.i) vW{gW=sW} = (^) M (i - £) 

The parameter range ]l,oo[ for p is called the supercritical regime, where 
there is an emerging giant cluster as N — > oo [4, 15]. Here, we say that a 
subgraph of Q( N > is a cluster, if it is a maximally connected subgraph of Q( N ) . 
By convention, we want to include isolated vertices as one-vertex clusters 
in this notion, too. In contrast to the supercritical regime, the subcritical 
regime p £]0,1[ has the property that the fraction of vertices which are 
either isolated or belong to tree clusters tends to one as ./V — > oo [4, 15]. 



LIFSHITZ TAILS FOR RANDOM GRAPHS 3 

Given any two different vertices i,j € {1, ...,N} of KL^ N \ i / j, let us 
denote the edge connecting i and j by the unordered pair We write 

ej^Q for the random variable which is one, if the edge is present in 
g( N ) . Otherwise, e\ N ) is zero. Hence, the e^'s are i.i.d. Bernoulli random 
variables with parameter p/N under the measure Pp^. 

Definition 2.1. The graph Laplacian A^ = A(Q^) of Erdos-Renyi 
random graphs Q( N *> is the random linear operator on with matrix ele- 
ments 

(2.2) a<;)=( £ ogjVejSIa-*,) 

for all i,j e {1, . . . ,N} in the canonical basis of C^. Here 5ij = 1 if i = j and 
zero otherwise denotes the Kronecker delta. 

Remark 2.2. (i) The diagonal matrix elements of specify the ran- 
dom vertex degrees in 

gW. if 

we denote the corresponding diagonal matrix 
by then (2.2) can be rewritten as 

(2.3) A W=#)-#), 

where A^ N ' is the adjacency matrix of the graph. 

(ii) The Laplacian A^) is nonnegative, as follows from its quadratic 
form 

N 

(2-4) (^A(%) = i^eg|^-^f 

for all tp G C^, where (•, •) stands for the standard scalar product in C^. 

(hi) We conclude from (ii) that the dimension of the kernel of A^-* is 
equal to the (random) number of clusters in — the corresponding eigen- 
vectors are constant within any cluster. 

(iv) Being a random self-adjoint and nonnegative N x iV-matrix, A^' 

possesses N nonnegative eigenvalues {^j}je{i,...,N} which are, of course, 
random variables. The normalized eigenvalue counting function 

(2.5) 4 N \E) := Ef )(iV- 1 #{j e{l,...,N}: \f> < E}) 

measures the average fraction of eigenvalues that do not exceed a given 
£?€l. Here, Ep^ denotes the mathematical expectation with respect to 
the probability measure 
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The quantity we are interested in is defined in the following: 

Lemma 2.3. Given any p £ ]0, oof, there exists a right- continuous distri- 
bution function a p : M — > [0, 1], which is called the integrated density of states 
of the Laplacian for Erdos-Renyi random graphs, such that 

(2.6) a p (E)= lim a p N \E) 

holds for all EeM, except for the at most countably many discontinuity 
points of o~ p . 

Remark 2.4. (i) The lemma is proven in Section 5, using the known 

fact [20] that the moments of a p N ^ converge as N — ► oo. An alternative 
approach to the proof of Lemma 2.3 via a resolvent-generating function was 
previously suggested in [20] . 

(ii) Remark 2.2(h) implies that a p {E) = for all E < 0. 

In this paper we focus on the subcritical regime p £]0, 1[, where there is 
no emerging giant cluster as N — > oo. Our main result is stated in the next 
theorem. 



Theorem 2.5. Let p £]0,1[. Then, a p has a Lifshitz tail at the lower 
edge of the spectrum, E = 0, with a Lifshitz exponent 1/2, that is, 

(0 7) ln\ln[a p (E)-a p (0)]\ _ 1 

{ ' ' eio In £7 2' 

Remark 2.6. (i) Theorem 2.5 follows from Lemmas 3.1 and 4.1 below, 
which provide upper and lower bounds for a p . Their proof is close in spirit 
to that in [23] for Laplacians on bond-percolation graphs. The bounds of 
Lemma 3.1 and Lemma 4.1 yield 

(2.8) p - 1 - lnp < - hm ^td^ZpM < 2 Vs(p - hip) 

for all E G]0,oo[ and all p g]0, 1[, which is a slightly stronger statement 
than (2.7). 

(ii) The main tool of the proof is a cluster decomposition of the block- 
diagonal Laplacian. The value 1/2 for the Lifshitz exponent relates to the 
fact that the asymptotic behavior of o~ p {E) as E \ is dominated by the 
smallest eigenvalues of the linear clusters in Erdos-Renyi random graphs. 
Indeed, the Cheeger-type lower bound of Lemma A.l for the smallest nonzero 
Laplacian eigenvalue of a cluster, which enters the upper bound for cr p (E) — 
crp(O), captures the correct size dependence for large linear clusters up to 
a constant. It also ensures that the smallest nonzero eigenvalue of a linear 
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cluster is among the smallest of nonzero eigenvalues of all clusters of the same 
size. Together with the exponential decay of the cluster-size distribution, this 
will yield the desired upper bound. On the other hand, the lower bound is 
obtained from retaining only the contribution of linear clusters to <J p (E) — 
cTp(O), which can be estimated in an elementary way. 
(hi) Remark 2.2(h) and (hi) imply that 

(2.9) 4 N) ^= * 

can be interpreted as the mean number density of clusters. Here, we intro- 
duced the orthogonal projection Pq on the kernel of and abbreviated 
the trace over N x iV-matrices by Tr. It is known [4] that the number of 
clusters grows linearly in N as TV — > oo, so E = has to be a discontinuity 
point of a p . The equality 

(2.10) <7 P (0)= lim 4 N \0) 

is therefore not guaranteed by Lemma 2.3, but we show in (3.9) below that 
it is true for every p g]0, 1[. It seems an open question to us whether (2.10) 
remains true for p G [1, oo[ — as is the fate of Theorem 2.5 for p G [1, oof. 

(iv) The papers [2, 20] provide recursion relations for the moments of 
the integrated-density-of-states measure as N — > oo for both the adjacency 
matrix and the Laplacian of Erdos-Renyi random graphs. The asymptotic 
behavior of the largest eigenvalue of the adjacency matrix was determined 
in [24], and [17] shows that the integrated density of states has a dense set 
of discontinuities; see also the numerical results in [2]. 

(v) More detailed spectral properties of Erdos-Renyi random graphs 
have been obtained in the theoretical-physics literature. Using the replica 
trick and other non rigorous arguments, it is argued in [6] that 

E10 £-1/2 

(2.11) 

= -[1 - p(l - Q p )} 1 / 2 ln[p(l - Q p )\ =: g(p) 

for all p G ]0, oof, where Q p is the biggest nonnegative solution of the equation 
Q = 1 — e -p< 2. If p G ]0, 1[, then Q p = and g(p) simplifies to —(1 — p) 1 / 2 lnp, 
which lies in between the bounds provided by (2.8). In [7, 14] the density of 
states of #) was examined by a combination of analytical and numerical 
methods for general E > 0. Their numerical results, however, were not con- 
clusive enough as to deduce the Lifshitz-tail behavior (2.7). The existence 
of emerging delocalized states in the giant cluster for p S> 1 was addressed 
in [5, 28]. 
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(vi) Weighted Erdos-Renyi random graphs also occur in some physical 
applications. The associated graph Laplacian is again given by (2.2), but 

now the probability distribution Pp^ of the edge random variables ef^.j is 
more general than Bernoulli. In many cases, it is required to have an atom 
at zero with weight 1 — p/N, corresponding to an absent edge, and a finite 

second moment Pp^{(ef^i) 2 } which is of the order iV" 1 . Random matrices 
of this type also fall under the name sparse random matrices. Some of their 
spectral properties were studied by, for example, [6, 18, 19, 28, 31]. 

3. Estimate from above. This section serves to establish an upper bound 
for the integrated density of states. The bound relies on a Cheeger-type 
inequality and the exponential decay of the cluster-size distribution in the 
subcritical regime. Both results are included in the Appendix. 

Lemma 3.1. There exists a constant C G]0,oo[ such that, for every p G 
]0, 1[, the integrated density of states satisfies the estimate 

e fip) 

(3.1) a p (E) - <7 P (0) < C expi-f^E- 1 / 2 } 

p 

for every E g]0, oo[, with the strictly positive decay parameter f{p) := p — 
1 — lnp. 

Proof. Fix p G]0,1[ and let E g]0,oo[ be a continuity point of a p . 
We introduce the right-continuous Heaviside unit-step function so that 
0(x) = 1 if x > and zero otherwise. Appealing to the spectral theorem and 
the functional calculus, we infer from (2.5) and (2.9) that 

a W( E ) - (rW( ) = Ef^N' 1 Tr[S(E - A^) - P (7V) ]} 

(3 ' 2) =Ef){[e(JS-AW)-P W] u }. 

To get the second equality in (3.2), we evaluated the trace in the canoncial 

basis and used enumeration invariance of Ep^ . 

The Laplacian has nonzero off-diagonal matrix elements between 

vertices in the same cluster only — and so has any function of A^'. Let C(l) 
denote the cluster of Q( N ' that contains vertex number 1. We define the 
associated graph Laplacian A(C(1)) as the random linear operator on 

whose matrix elements [A(C(l))]y coincide with A^^ for i,j G C(l), but are 
zero otherwise. Likewise, Po(C(l)) stands for the orthogonal projector in 
on the kernel of A(C(1)). Thus, introducing the characteristic function Xfi 
of the event that vertex number 1 is not isolated, we obtain 



af)(E)-af)(O)=Ef){xa o [0(^-A(C(l)))-P o (C(l))] 11 } 
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(3.3) <Ef){xn o e(£-£ mi n(C(l)))} 

<Ef){ XQo e(£;-|c(i)|- 2 )}. 

The first inequality in (3.3) follows from the spectral theorem with E m i n (C(l)) 
denoting the smallest nonzero eigenvalue of A(C(1)). The second inequality 
in (3.3) uses the Cheeger-type estimate of Lemma A.l in the Appendix, and 
|C(1)| counts the number of vertices in the cluster C(l). The expression in 
the last line of (3.3) is equal to 

oo m(E)-l 

(3.4) ^{\C(l)\=n} = l- £ Pf){|C(l)|=n}, 

n=m(E) n=l 

where m(E) :=max{2, L-^' _1//2 J} an d [x\ stands for the biggest integer, not 
exceeding iGR, Since E was chosen to be a continuity point of a p , we 
deduce with the help of Lemma 2.3 that 

m(E)-l oo 

(3.5) a p (E) -limMaW(0)<l- V nr n {p)= T nr n { V ). 

n=l n=m(E) 

Here we introduced the cluster-size distribution 

(3.6) th(p) := n- 1 lim Pf ) {|C(l)| = n} 

N— >oo r 

of Erdos-Renyi random graphs, whose existence (A. 7) and normalization 
(A. 8) is summarized in Lemma A. 2 in the Appendix. Taking the limit E J. 
in (3.5) along a sequence of continuity points and appealing to the right- 
continuity of a p , we infer 

(3.7) lmiinfaf)(0)>a p (0). 

On the other hand, the monotonicity and right-continuity of a p imply, for 

all p £ ]0, oof, 

(3.8) lim sup (7^(0) < lim lim sup ^(E) = \\ma p {E) = a p (0), 

where the limit E j is again taken along a sequence of continuity points. 
From (3.7) and (3.8), we conclude the existence of the limit 

(3.9) lim <7W(0) = ffp (0) 

N— >oo 

in the subcritical regime pe]0,l[. Hence, (3.5) and the exponential decay 
(A. 9) of the cluster-size distribution lead to 

1 oo 

3/2 -n/(p) 



a p (E)-a p (0)<^=- J2 n ~ V2e 
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-m(E)fip) oo 



(3.10) < 6 V [n + mCE)]- 3/2 e- n/ W 

V2tt P n=Q 



^ exp{-/(p)[E-V 2 -l]} °° _ 3/2 

- ./ST.. 



n=2 



Finally, the estimate (3.10) extends to all E E ]0, oo[ by right-continuity. This 
completes the proof of Lemma 3.1. □ 

4. Estimate from below. In this section we derive a lower bound for 
the integrated density of states by retaining only contributions from linear 
clusters. 

Lemma 4.1. Let p £ ]0, 1[ and define F{p) :=p — Inp > 1. Then the es- 
timate 

(4.1) a p (E) -o-p(O) > — — exp{-2 v / 3F(p) J E- 1 / 2 } 

ZioWs /or every E £ ]0, oof. 

Remark 4.2. The decay parameter .F(p) is related to that in Lemma 3.1 
by F(p) = f(p) + 1. 

Proof of Lemma 4.1. Fix p &]0, 1[. The right-continuity of a p and 

the monotonicity of a p N ^ imply a p (E) > limsup N _ too a p N \E) for every E 1 € 
]0,oo[. Together with (3.7) and (3.2), this yields 



(4.2) <J P (E) - Op(0) > limsupE^lAf- 1 Tr[9(£ - A^) - P W ]}. 

AT— ►oo 

The decomposition of the random graph Q( N > into its random clusters 
e{i,...,K} provides us with the relation 

Tr[@(E - AW) - P (7V) ] = J2 ^ME - A(cf } )) - P (cf°)] 

fe=i 

(4.3) 

>j;9(|4"l|-2)9(E-B„ lin (C«)). 

fc = l 

The inequality in (4.3) relies on the spectral theorem. Note that one-vertex 
clusters do not contribute to the right-hand side of the first line in (4.3). 

For n £ {2, . . . ,N}, let X£„(£-£ ) be the indicator function of the event 

that the cluster c£ is a linear chain with n > 2 vertices, that is, that it is 
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a connected graph having n — 2 vertices with degree 2 and 2 vertices with 
degree 1. Using (4.2) and (4.3), we then obtain the first inequality of the 
chain 

a p (E) - Op(0) 

r K oo 

> hmsupEf ) iV- 1 E £ X£n ( C f >)9(i5 - ^(cf } )) 



v- 



(4.4) 



fe=ln=2 



> 



limsup £ 0(£ - 12/n 2 )Ef ) iV" 1 £ X£„ (Q 



v 



k=l 



> — - — - limsup ¥^{C(1) is linear and has M(E) vertices}. 

M(E) N^oo 

To derive the second inequality in (4.4), we used the upper bound 12/n 2 for 
the smallest nonzero Laplacian eigenvalue of a linear chain with n vertices, 
see, for example, Lemma 2.6(i) in [23]. For the last inequality in (4.4), we 
introduced M(E) := [(12/E) 1 ' 2 \ + 1, the smallest integer strictly greater 
than (12/ E) 1 / 2 , dropped all terms in the n-sum, except the one with n = 
M(E), and observed 

Ex £ Jcf))) = ^Ef){f: X£ri (C( J ))) 



(4.5) 




~Ef){ X£n (C(l))}. 
n y 



Equation (4.5) involves C(j), the cluster of containing vertex j, and 

it exploits enumeration invariance of Fp . Now, elementary combinatorics 
shows for any m G 2, . . . , N that 

¥^{C(1) is linear and has m vertices} 

(4.6) 

/ N _ 1 \ m! / p\ m ^ / p\ (JV-3)(m-2)+2(AT-2) 
= \m-l) ~2\n) \~n) 
Here the first factor corresponds to the choice of m — 1 vertices different 
from the already fixed vertex number one. The second factor corresponds 
to ordering these m vertices in a chain. The third factor accounts for the 
probability to join these m vertices by m — 1 edges and the last one assures 
that there are no other edges joining the m vertices to the remaining N — m 
vertices. Hence, the limit N — > oo of (4.6) exists and is given by 

lim ¥^{C(1) is linear and has m vertices} 

TV-^oo 

(4.7) 

= ™ m-l e ~pm Um (iV-1)! = m m _i e _ pm 

2 F N^ooN m ~ 1 (N -m)\ 2 1 
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Inserting this into (4.4), we arrive at 

(4.8) a p (E)-a p (0) > ±- exp[-(p - hxp)M(E)], 

which implies the lemma. □ 

5. Existence of the integrated density of states. 

Proof of Lemma 2.3. Let p e]0, oo[. Theorem 2 in [20] establishes the 
existence and finiteness of the limits 

(5.1) M fe A := lim / da { D N \E)E k = lim ^{N" 1 Tr[(A^) fc ]}, 

AT->ooJ[0,oo[ iV->oo 

k £ No, of all moments of a p N ^ as N — > oo. Being the limit of a sequence of 
Stieltjes moments, {M ^}k^n is itself a sequence of Stieltjes moments as- 
sociated to some, not necessarily unique, distribution function a p on [0,oo[. 
This follows from Theorem 1.1 in [32], see also the statements in Chap- 
ter 1.2(b) there. We extend a p to M. by setting it to zero on ]— oo, 0[. We will 
prove the bound 

(5.2) M 2 A fc < (c p k) 2k 

for all k £ N with some ^-independent constant c p G]0,oo[. This, in turn, 
guarantees the Carleman condition Y^k=i(M^)~ l ^ 2k ^ = +oo and, by Theo- 
rem 1.10 in [32], the uniqueness of the Hamburger (and, hence, the Stieltjes) 
moment problem. Knowing the uniqueness of a p , the lemma then follows 
from, for example, Theorem 4.5.5 in [11]. 

To prove (5.2), we use (2.3) and Holder's inequality for the Schatten trace 
norms \\B\\ q := {Tr[(B*By/ 2 }} 1 / q , q > 1, of complex N x iV-matrices, where 
B* stands for the adjoint of B. For every natural number N > p, this yields 
the bound 



Tr[(AW) 2fc ]<^( 2 ^)||^)||^pW||^ 



K=0 

2k 



(5-3) =(\\D^\\ 2k + \\A^\\ 



2k) 

<2»- 1 (||DW||2 + ||AW|a 



Hence, we get 

(5-4) M 2 i<2 2fe - 1 (M 2 D fc + M 2 l), 

where, thanks to ergodicity, 

oo n 

(5.5) MR := lim E^AT 1 Tr[(D^) 2k ]} = e~ p V V —n 2k 

N-,00 ^ nl 
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is nothing but the 2fcth moment of the Poissonian [4] vertex-degree distri- 
bution and := limn^E^ {N' 1 Tr[(A^) 2k }} is the 2fcth moment of 
the adjacency matrix. So, (5.2) is implied by (5.4), provided we show 

(5.6) max{M&,M&} < (c p k/2) 2k 

for all k G N. Concerning M^, this follows from applying the elementary 
inequality a b < b b + b a to (5.5). This inequality holds for a G ]0, oo[ and b G 
[e, oof. Its validity is obvious for a<b, while, for a > b > e, it can be deduced 
from Jensen's inequality. The desired bound for is established in the 
proof of Proposition 1 in [17], using a result of [21]; see also [2]. Hence, (5.2) 
is proven, and so is the lemma. □ 

APPENDIX 

For completeness and convenience of the reader, we state and prove two 
auxiliary results in this Appendix, which were needed in the proof of the 
upper bound in Lemma 3.1. The first result concerns a weakened version 
of a Cheeger-type inequality, which does not involve the graph's maximum 
vertex degree. 

Lemma A.l. Let C be a connected finite graph with \C\ > 2 vertices. 
Then the smallest nonzero Laplacian eigenvalue E m [ n (C) is bounded from 
below according to 

1 



(A.l) E min (C) > u 



Proof. The proof is inspired by that of Lemma 1.9 in [10]. It does not 
involve the maximum vertex degree, though. 

Set n := \C\ and let V := {1, . . . , n} be the vertex set and £ the edge set of 
C. Elements of £ are denoted by unordered pairs [i,j] of the vertices i,j G V 
they join. The minmax-principle and (2.4) imply that 

(A.2) E min (C)= inf E[ ^7 j)2 > 

where the infimum is taken over M. n only (instead of C n ), because all eigen- 
vectors of A(C) can be chosen to be real. The other constraint expresses the 
fact that the nondegenerate zero eigenvalue of the connected graph corre- 
sponds to an eigenvector with constant components. 

Now, for any given (p G R n , obeying the orthogonality constraint J2i^v V 9 * = 
0, let u G V be such that \<p u \ = maxjgy \ ipi\. Due to the constraint, there ex- 
ists v G V such that 

(A.3) (p u <p v < 0. 
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Let C £ be the shortest path in C connecting the vertices u and v. Then 
we have 

(A.4) E min (C) > mf — — . 

The triangle and the Cauchy-Schwarz inequality supply us with the estimate 

(A. 5) \<p u -<Pv\< \<Pi-<Pj\^\ E (Vi-^jff l^| 1/2 > 

[i,j]ev v <[i,j]er v ) 

where \P<p\ stands for the number of edges in Vu>- Inserting (A. 5) into (A.4) 
and noting \V V \ < \C\, we arrive at 

(A.6) £„,„(C)>-L i„ f rf + rf^ , 

\C\%m™:J2 iev <Pi=o ft 

The claim now follows from (A. 3). □ 

The second auxiliary result summarizes the existence, normalization and 
decay of the cluster-size distribution of Erdds-Renyi random graphs in the 
subcritical regime. 

Lemma A. 2. Assume the subcritical regime p e]0, 1[ and let C(l) be the 
maximally connected subgraph of QW) containing vertex number one. Then 
the cluster-size distribution 

(A.7) rjp) := n~ x lim pW{|C(l)| =n} = i-n n "V~ 1 e~ np 

exists for every n E N and equals the mean number density of tree clusters 
with n vertices. It is normalized according to 

oo 

(A.8) £ nr n (p) = 1 

n=l 

and has an exponentially small tail 

( A -9) T n (p)<-^^j- 2 e- n ^\ 

with the decay parameter f(p):=p — l — hip > 0. 

Proof. The lemma follows from collecting some well-known properties 
of Erdds-Renyi random graphs in [4]. Fix n G N and let C(l) S T, respec- 
tively C(l) G T n , denote the event that is a tree cluster, respectively, a 
tree cluster with n vertices. Then we have 

Pf){|C(l)|=n} = Pf){C(l)eT n } 

(A.10) 

+ pW{|C(l)|=nandC(l)gT}. 
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Introducing \A, the characteristic function of an event A, we get an upper 
bound for the last probability in (A. 10): 

p(*>{|C(l)|=nandC(l)£T} 

(A.ll) <PW{C(1)£T} 

1 N 

= 1 - Pf ){C(1) E T} = 1 - - ^Ef ){ Xr (C(i))}. 

The sum in the last line of (A.ll) represents the mean number of vertices 
on tree clusters in Q^ N \ Hence, Theorem 5.7(h) in [4] implies that 

(A.12) lim PW{|C(1)| = n and C(l) i T} = 

for all pG]0,l[. On the other hand, recalling the notation for enumerating 
clusters above (4.3), we deduce from the equality 



1 1 N 

-Pf ){C(1) E T n } = — ^Ef ){xr n (C(j))} 



(A,3) f 

= ^E<Wcf°)} 
fc=i 

that n" 1 pJ ) Ar) {C(l) E T n } equals the mean number density of tree clusters 
with n vertices in Q^ N >. Accordingly, (5.1) in [4] yields 



^w.)^.)4(I)»'1i 

(A.14) 



n-1 



x n(N-n)+[ _ )-n+l 



X 1- — 

This expression has a limit as N — > oo, which is given by 
(A.15) r n (p) := - lim Pf ){C(1) E T n } = 

Taken together, (A. 10), (A.12) and (A.15) establish the first assertion (A. 7) 
of the lemma. The second assertion, the normalization (A. 8) follows from (5.6) 
in [4] and the definition in the equation above (5.5) in [4]. Finally, to prove 
the decay (A. 9) of T n (p), we apply the Stirling inequality n\ > (n/e) n \/2vrn x 
exp{l/(12n + 1)}, see, for example, (1.4) in [4], to (A.15). □ 
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